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Abstract 

Dichotomous noise appears in a wide variety of physical and mathematical models. It has escaped 
attention that the standard results for the long time properties cannot be applied when unstable 
fixed points are crossed in the asymptotic regime. We show how calculations have to be modified 
to deal with these cases and present as a first application full analytic results for hypersensitive 
transport. 

PACS numbers: 05.40.0a,02.50.-r 



While the Wiener process together with its "time derivative," the Gaussian white noise, is 
certainly the method of choice to describe Brownian motion, the motion induced by another 
fundamental stochastic process, namely, the dichotomous Markov process (see, e.g. 01); 
has its own virtues and interest. Systems driven by dichotomous noise can often be de- 
scribed analytically, including the Gaussian white noise case as a specific limit, and allow 
the analytic investigation of the effects of the finite correlation time of the noise, notably in 
noise induced transitions, noise induced phase transitions, stochastic resonance, and ratch- 
ets. Dichotomous noise is often a good representation of the actual physical situation, e.g., 
thermal transitions between two configurations or states, and can easily be implemented as 
an external noise, with the additional advantage that the support of this noise is finite. 

A widely studied generic stochastic equation that describes the temporal evolution of a 
single scalar variable x{t) is 



where the dot stands for the time derivative, ^(t) is a symmetric dichotomous Markov process 
that takes on the values ±1 with transition rate k, v{x) is a given velocity profile, and F is 
a constant external force. One can of course embellish this description in a variety of ways 
such as, e.g., by allowing for a state- and/or time-dependent external force, but here we 
adhere to this simple form. Existing results include the steady state distribution P, ^ and 
first passage time moments, see e.g. 0. When Eq. (|ID is defined on an interval with periodic 
boundary conditions, one can evaluate the stationary probability fiux and from it the average 
asymptotic drift velocity or diffusion coefficient, a problem that has recently received a great 
deal of attention in the context of Josephson junctions and Brownian motors 0. Although 
these results are often claimed to be completely general, our study shows that this is not 
the case. Indeed, to our knowledge, with the exception of [0, all the existing results are 
limited to motion that asymptotically does not cross unstable fixed points of the dynamics. 
The main purpose of our work is to point out where the existing results break down, and to 
present the procedure to obtain a fully general solution for the asymptotic average velocity, 
including, as a first direct application, the problem of hypersensitive response [§]. 

Consider, then, the prototypical stochastic differential equation (|l|). We take the velocity 
profile v{x) to be periodic, v{x) = v{x + L), with zero average, v{x) dx = 0. A schematic 
representation of the two configurations assumed by the "net potentials" ^V{x) — Fx asso- 
ciated with right hand side (with v{x) = —dV{x)/dx) is shown in Fig. |l]. The fixed points 
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FIG. 1: The net potentials =fF(x) — Fx. 



of the dynamics are the points at which ±f (a;) + F vanish, i.e., the local extrema of the 
net potentials. The stochastic dynamics (|l]) can equivalently be described by the master 
equation for the probability densities P+{x, t) and P^{x, t) for being at x at time t if ^ = +1 
and —1, respectively, with x G [0,L] and periodic boundary conditions: 



To find analytic expressions for the asymptotic steady state (i.e., time-independent) prob- 
abilities and mean velocity (x), it is more convenient to introduce the sum and difference 
probability densities P{x, t) = P+{x, t) + t) and p{x, t) = P+{x, t) — P-{x, t). Summa- 
tion of the equations in then immediately leads to the conclusion that in the steady state 
the probability flux J associated with P{x), namely, J = FP{x) + v{x)p{x), is a constant 
whose value is to be determined. It also leads to a direct relation between the mean velocity 
in the stationary state and the flux, {x) = {[f (x) + F] P+{x) + [— f (x) + F] P_(x)} dx = 



where v'{x) = dv{x)/dx. The solution to this equation, together with the constant flux con- 



dP±{x,t) 
dt 



Tg^{[v{x)±F]P±{x,t)} 
-k[Pi{x,t)-P^{x,t)]. 
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dition and the normalization condition for P{x), P{x)dx = 1, can be used to determine 
P{x), the flux J and the mean velocity (x). 



The crux of the problem now resides in finding the solution to Eq. (j^). This solution 
is straightforward when [F^ — v'^{x)] has no zeroes, that is, when the net potentials have 
no extrema within the interval (0,L). In this case the standard method of variation of 
parameters leads to the familiar solution 

Pix) = -JpT^:;^^ [CGix, 0) + Hix, 0; x)] , (4) 

where C is a constant of integration that arises from the general solution to the homogeneous 
part of Eq. (H), the second contribution is the particular solution, and we have defined the 
functions 

H{z,y;x) = f sgn[F'^ - v'^{x')\v' {x')G{x,x')dx\ 

J V 



G{z,y) = exp 
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The usual procedure to determine C is to require periodicity of p{x), recalhng that v{x) is 
periodic. One finally obtains p{x) = [J — FP{x)]/ Jv{x), and 



F [ [1-G{0,L)][F^ -v^{x)] 
X / dx'v\x')G{x,x') \ . (6) 
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The normalization of P{x) determines the value of the fiux J and leads to the following 
result for the mean velocity at the steady state: 



^[1-^(0,^)] 7o F^~v^{x) 

x+L ^ -1 

dx'v'{x')G{x,x')\ . (7) 

The standard results shown above are applicable not only in the absence of fixed points, 
but also when the asymptotic behavior is governed by stable fixed points. In this latter 
case the dynamics settles into an alternating motion between these points, so that they 
delimit the interval in which the steady state probability is non-zero 0. The associated 
normalizable divergences at the fixed points represent regions where the probability density 
for finding the system is high. 

The situation is entirely different, both physically and mathematically, when the system 
can cross unstable fixed points within the interval (0, L) in the long time limit. A simple 



illustration is provided by the example v{x) = sin x. In the absence of an external force, the 
dynamics is restricted to an interval [kn, (/c + l)7r] {k integer). Even though the application 
of an external forcing \F\ < 1 cannot induce "escape" from this interval in either of the 
separate dynamics x = sin x + F and x = — sin x + F, running solutions with finite average 
velocity appear when the dynamics switches back and forth between the two (see Fig. p. 
The explicit calculation of this velocity is one of our main goals. Clearly, the solution (j^) is 
no longer correct because it contains non-integrable singularities (see below) at the unstable 
fixed points where the probability of finding the system is expected to be low, not high. 

For simplicity we restrict our presentation to velocity profiles v{x) that are continuously 
decreasing functions of x on [0, L/2] and symmetric about L/2, v{x + L/2) = —v{x). This 
implies that P{x + L/2) = P{x) and p{x + L/2) = —p{x), so that we can limit our analysis 
to the interval [0, L/2]. In this "minimal scenario", the equation F"^ — v'^{x) = has only 
two solutions in the interval [0, L/2], namely, Xi, corresponding to an unstable fixed point 
in the ^ = — 1 dynamics [F = v{xi)], and X2, a stable fixed point in the ^ = +1 dynamics 
[F = —v{x2)], with X2 > xi. The steady state results leading to Eq. (^ still apply, but the 
solution to Eq. is more delicate than the "blind" integration that yields Eq. (|^). Indeed, 
the coefficient of the first derivative is zero at the fixed points, which now lie entirely within 
the support of the probability distribution. Thus, the equation becomes singular. 

The method of variation of parameters for an equation of the type (^ leads to a solution 
which is a sum of the general solution of the homogeneous equation and a particular solu- 
tion of the inhomogeneous equation, as in Eq. (|^). The subtlety lies in the determination 
of the constant of integration C, which in the previous case was fixed simply by imposing 
periodic boundary conditions. In the vicinity of the stable fixed point X2 this straight- 
forward procedure leads to the dependence P{x) ~ \x — X2\''^^"'^^^'^^~^, which is continuous 
when k/\v'{x2)\ > 1 and divergent but integrable for k/\v'{x2)\ < 1. This result causes no 
mathematical difficulty and is consistent with the physical intuition that probability near 
a stable fixed point should indeed build up, especially when the switching rate is low. At 
the unstable fixed point xi, however, this procedure leads to an apparent non-integrable 
divergence, P{x) ~ |x — xil^'^/''"''-^^^'"^, which is clearly unphysical and mathematically im- 
proper in view of the requirement of normalization. The fallacy lies in the assumption that 
a single constant C is valid throughout the region (0,L). In fact, the solution (^ is valid 
in the separate intervals [0,a;i), {xi,X2), and {x2,L/2], but not necessarily with the same 
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constant of integration in all of them. Indeed, there is exactly one choice of this constant 
valid for both [0,a;i) and {xi,X2) such that the divergence at xi is removed, and another 
choice valid in the interval {x2,L/2] that ensures required continuity and periodicity. In 
other words, even though the general solution of the homogeneous equation always has a 
divergence, there exists exactly one solution to the full inhomogeneous equation that has 
no divergence and is actually completely smooth at xi. This solution is given by Eq. with 
the choice C = —H{xi, 0; 0) in both intervals [0, xi) and (a;i, X2) 0. This choice insures that 
the coefficient of the divergent term vanishes at x = Xi. We conclude that for x G [0,^2), 

J 



v[x , 



(8) 



Note that P{x) is now continuous at Xi, and that lima-^^.^ P{x) = lim^./^a,^ P{x) = 
JF-^ {1 - [2{k/\v'{xi) \ + 1)]-^} is indeed finite. 

For X G {x2,L/2], the result (|[) for p{x) applies again, but now the constant C is 
determined by imposing the continuity of p{x) at x = L/2. One finds: 

+G(0,L/2)i7(xi,0;x)]}. (9) 

At the stable fixed point X2, P{x) has the behavior described earlier, i.e., it is continuous 
for k/\v'2{x2)\ > 1 and divergent but integrable for k/\v'{x2)\ < 1. 

The values of the fiux J and of the average velocity follow from these results by imposing 
the normalization of P{x): 

2 f^^"^ v(r) 
+ 7 / dx-^^^ [H{x,L/2;x) 

+G(0,L/2)if(x,0;x)]}"\ (10) 

This is our main new result. Note that the above procedure can be repeated straightfor- 
wardly but tediously for more complicated cases involving several stable and unstable fixed 
points. 

To illustrate our findings with explicit results, we turn to a particular case of a piecewise 
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FIG. 2: Probability density P{x) vs x/L for the parameter values / = 0.5, T = 0.4, and a = 1.0. 
Histogram: simulation results. Curve: exact theory. 



linear profile: 



ViX) 



vo, for X e [0, L/2 -21) 

vo (^^ - 1 - jY ioT X e [L/2 - 21, L/2) 



-v{x - L/2) 



for X e [L/2, L) 



(11) 



with I < L/A and, of course, the periodicity condition v{x + L) = v{x). It is convenient to 
introduce the following dimensionless variables: 



/ = F/vo, a = Ik/vo, r = 4//L. 



(12) 



In this case the function H 


[x,0;0)/vo 


= T{x 


for X e [L/2 -21, L/2), 








T{x) 
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(13) 



Explicit exact results for the probability densities and for the resultant average velocity for 
all values of / are available and will be detailed elsewhere P]. Here we exhibit only some 
of these results for the new case < / < 1. Figure |^ depicts a typical probability density 
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FIG. 3: Mean velocity as a function of the applied force for various values of a. Note the 
hypersensitive response. 

P[x) MS x/L that clearly shows the agreement between the exact theoretical results and 
simulations. For the average velocity we find 



The above integrals can be evaluated explicitly for specific values of a, in particular, 



slow switching rate, that is, in the adiabatic regime a — 0, one recovers the region of 
hypersensitive response discussed in namely, (i) ~ 2vQa/T. The physics of this result is 
explained as follows: The stable and unstable fixed points of the dynamics, which coincide 
in the absence of forcing, are shifted apart by an amount of order of IF/vq for a small force. 
Upon each switch of the dichotomous process and for sufficiently slow switching rate, the 
particle will glide down to the next stable fixed point, crossing the location of the unstable 
fixed point of the alternate dynamics, cf. Fig. ^ As the average time between switches is 




(14) 



for a 
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k~^, and the distance covered is L/2, the mean velocity is just (x) = Lk/2, independently 
of F and fo- The typical time r for a particle to escape the region around the fixed point is 
determined by the relation (IF/vq) exp(for//) ^ /. The crucial observation for hypersensitive 
response is that the necessary condition, r << k^^ or a << {—Inf)^^, can be satisfied by 
very small forces for moderately small a because of the logarithmic dependence on /, see 
Fig. ^. In this figure one also observes the region of "normal" (i.e., linear) response for higher 
forcing or frequency and, more relevant to our preoccupation here, the strongly nonlinear 
dependence at very low forcing. In fact in the limit / — 0, one finds from Eq. (14) that 
{x)/vo ~ [2ar(/n/)^]~^. In other words, hypersensitive response is very pronounced in this 
region, with the velocity picking up with an infinite derivative at / = 0. 

In conclusion, the procedure presented here has resolved all technical problems related 
to steady state dichotomous dynamics, making possible the analytic description of cases 
involving the crossing of fixed points in the asymptotic regime. 
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